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ATTENUATION OF SINUSOIDAL DISTURBANCES IN NONVISCOUS LIQUID 
FLOWING IN A LONG LINE WITH DISTRIBUTED LEAKAGE 
by Robert J. Blade and Carl M. Holland 
Lewis Research Center 

SUMMARY 

Sinusoidal disturbances from 0. 5 to 82 cps superimposed on the flow of a nonviscous 
liquid were attenuated to design values in a 66 -foot-long (20. 1 m) line. The disturbances 
were attenuated by placing shunt assemblies that simulate distributed leakage along the 
line. Experimental values of the attenuation constant were determined and compared with 
the design values that were obtained from analogous electrical transmission line theory. 
The results showed that sinusoidal perturbations at a given frequency can be attenuated 
to a design value by using shunt assemblies and that theory and experiment agreed well. 


INTRODUCTION 

Because acoustic pressure and flow perturbations are often of importance in the op- 
eration of various fluid systems, the Lewis Research Center has been involved in a gen- 
eral experimental research program to determine the effects of acoustic disturbances 
in flowing -liquid systems. Transfer functions for acoustic disturbances superimposed 
on the mean flow of liquid propellant in long lines of several geometric networks are 
reported in references 1 to 3. A study of the attenuation constant associated with sinus- 
oidal perturbations superimposed on the mean laminar flow of a viscous liquid in a long 
line is described in reference 4. 

Pressure and flow perturbations in long liquid-filled lines are often a contributing 
element toward instabilities in liquid rockets and hydraulic control systems. For exam- 
ple, longitudinal instabilities in liquid rockets have resulted from a closed-loop coupling 
of disturbances in the propellant feed lines, engines, and structure. A convenient meth- 
od of stabilizing these systems is to attenuate sufficiently the pressure and flow pertur- 
bations in these lines. 

This investigation was conducted to determine the efficacy of attenuating to design 


values sinusoidal disturbances superimposed on the mean flow of a nonviscous liquid 
(the type of fluid often used in practice) in a long line by means of shunt assemblies 
placed along the test line. These assemblies were closely spaced relative to disturbance 
wavelength. Experimental values of the attenuation constant were determined and com- 
pared with the values obtained from the analogous electrical transmission line equations 
for the zero distributed series disturbance -resistance and nonzero disturbance - 
resistance and nonzero distributed leakage. 

Attenuation studies were conducted with JP-4 fuel flowing in a 66 -foot-long (20. 1 m) 
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stainless -steel line (1 in. o. d. or 2. 54X10” m). The line was firmly supported at both 
ends and contained 33 equally spaced shunt assemblies. Two types of shunt assemblies 
were tested. One type was essentially a porous -metal disk. The other type was essen- 
tially 21 small-diameter tubes (lossy tubes) arranged in parallel. Flow was modulated 
at frequencies from 0. 5 to 82 cps by means of an electrohydraulic servothrottle. Mean 

line pressures ranged from 75 to 155 pounds per square inch gage (5. 17x10 to 
6 2 

1. 07x10 N/m gage), and mean line flow ranged from 7 to 9. 4 feet per second (2. 1 to 

2. 9 m/sec). 


APPARATUS 
Flow System 

The essential parts of the open-loop pumped-return flow system used in the experi- 
ment are shown in figure 1(a). The liquid (JP-4 fuel) was forced through the test line by 
a gear pump, and the mean flow rate was measured by means of a rotameter. A heat 
exchanger in the fluid supply tank maintained the fluid at 79°±3° F for the porous-metal 
shunt element studies and at 85°±3° F for the lossy-tube shunt-element studies. Hydrau- 
lic accumulators were placed between the pump and the test line to provide steady supply 
pressure. The discharge from the test line was submerged in a constant-height vented 
tank. Fluid was returned to the supply tank by intermittent operation of the return pump. 


Test Line 

The test line consisted of a 66-foot-long (20. 1 m) stainless -steel tube with a 1.00- 
inch (2. 54x10 m) outside diameter and a 0. 065-inch-thick (1. 6x10” m) wall, which 
was fitted with 33 equally spaced shunt assemblies as shown in figure 1(b). The test line 
rested horizontally on transverse wires spaced at 2 -foot (0. 6 m) intervals. The upstream 
and downstream ends of the tube were firmly supported in order to avoid line vibration. 
The downstream end was terminated in an orifice plate containing 34 holes of 0. 040-inch 

2 


■■■■■■■■III Ill HI 



r-Electrohydrauiic 




_ 3 

(1. 0x10 m) diameter. The range of run-averaged dimensionless downstream admit- 
tance for the runs was from 0. 95 to 1. 4. (Dimensionless downstream admittance is the 
ratio of downstream admittance to characteristic admittance of the line. ) 


Shunt Assemblies 

Schematic drawings of the two types of shunt assemblies used are shown in figure 2. 

O 4 o 

For both types of assemblies, an accumulator (total volume, 10 in. or 1. 64xl0~* m°) 
was used as a sink and source of liquid. Provision was made for bleeding to ensure that 
the fluid flowing through the shunt element was JP-4 fuel. 

The porous -metal shunt element was a commercial bronze, porous -metal disk 
(Morain Products Grade 4 Porex) 0.0156 foot (4. 75X10" 3 m) in diameter (the same as 
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(a) Porous metal. 


H.d., 5/8 in. (1.59xl0" 2 m) 

v r Accumulator 10 in? (1. 64xl0“ 4 m 5 ) 

\ ' total volume 

\ \ rBladder 


(b) Lossy tube. 

Figure 2 . - Typical shunt assemblies. 

Q 

the flow diameter) and 0. 0208 foot (6. 34x10” m) thick. It was mounted in a hole drilled 
through a stainless -steel block. 

The lossy-tube shunt element was essentially 21 parallel Inconel tubes of 0. 125-foot 
(3. 81x10”^ m) length and 6. 12x10”^ foot (1. 87x10”^ m) inner diameter mounted in a . 

6. 55x10” ^ -foot-diameter (2. 00x10”^ m) hole drilled through a stainless-steel block. The 
volume between these lossy tubes was filled with solder at both ends of the stainless -steel 
block. 

Flow Disturbance Generator 

Sinusoidal perturbations of flow and pressure (from 0. 5 to 82 cps) were induced in 
the system by means of an electrohydraulic servothrottle located at the upstream end of 
the test line. The throttle was oscillated sinusoidally (by a master oscillator) about a 
partly open mean position in response to an alternating voltage. The amplitude of the 
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sinusoidal pressure perturbations was of the order of 10 pounds per square inch 
(7xl0 4 N/m 2 ). 


Instrumentation 

Commercial flush-diaphragm pressure transducers were used to measure pressure 
perturbations at stations A, B (chosen for convenience), and C (fig. 1(b)). The trans- 
ducers with their associated electronic equipment were statically calibrated to determine 
their gain factors. The inphase and quadrature components (relative to the master 
servothrottle oscillator) of the pressure perturbations were determined by means of a 
commercial transfer-function analyzer, as was done for previous studies (refs. 1, 2, 
or 4). 


PROCEDURE 

Test 

Each accumulator associated with the shunt assemblies was charged to about 80 per- 
cent of the desired mean line pressure. Then the flow system was operated until the 
flow rate, temperature, and line pressure had stabilized. The mean flow rate, the pres- 
sures at stations A and C, and the liquid temperature adjacent to the orifice in the re- 
ceiver tank were recorded. The servothrottle was operated at a series of frequencies. 
The amplitude of the throttle -area variation was maintained constant over the entire 
range of frequencies for each run and was small relative to the mean open area in order 
to avoid nonlinear effects. At each frequency the values of the inphase and quadrature 
components of the pressure perturbations at stations A, B, and C (fig. 1(b)) were re- 
corded. The mean pressure drop across the downstream orifice was constant during 
each test run but was changed from run to run (by changing the mean line pressure) to 
vary the orifice impedance. 


Data Analysis 

The following quantities were computed from the experimental data for each fre- 
quency: 

(1) Attenuation constant, o', Np/ft 

(2) Phase constant, /3, rad/ft 

(3) Complex dimensionless admittance at station C, z/Zq (fig. 1(b)) 


r 
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(All symbols are defined in appendix A. ) 

The values of these quantities at each frequency were obtained from the following 
(complex) equations: 


p 

— = cosh(o! + j£)f A + -5_ sinh(a + j/S)f A (1) 

P C Z C 


— = cosh(a + j0)f B + — sinh(a + B (2) 

p c z c 

This pair of equations was obtained from the well-known pressure transfer function for 
damped sinusoidal one -dimensional waves in a uniform medium. (This transfer function 
is considered herein to be the defining relation for a. ) Equations (1) and (2) were 
solved numerically on a digital computer by a procedure described in reference 4. The 
necessary input data were (1) the six measured values of the inphase and quadrature 
components of the pressure disturbance amplitudes at stations A, B, and C (fig. 1(b)), 
that is, P A’ P B’ and P c ; and (2) the distance in feet between stations A and C, and 
between stations B and C, SL g. 


RESULTS AND DISCUSSION 


The experimental results were compared with analytical values obtained from the 
analogous electrical transmission line equation for nonzero distributed leakage and zero 
distributed series disturbance -resistance. As shown in appendix B, the analytical rela- 
tions for the attenuation constant and the phase constant, respectively, are given by the 
positive roots of 


“i = 




(3) 


and 


fi 


o> 2 K d L d + 



( 4 ) 


For the porous-metal shunt-element studies, the effective distributed leakage has 
the constant value given by 
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( 5 ) 


n _ 1 n 
° " R s 1 A 


and Kq and have the constant values (see appendix C) for zero shunt admittance. 
For the lossy-tube shunt-element studies, G D and K d are frequency dependent as 
given by 


G D " 


n 


r t 

_2 2 t 2 
Rrp + U> Lrp 


- — 


D ■ D 


£. d 2 ,, 2 t 2 

A R + (0 L 


( 6 ) 

( 7 ) 


and Lq again has its zero-leakage constant value. (Appendix C contains the numerical 
values of the parameters. ) 

The experimental values obtained for both the porous-metal and the lossy-tube shunt- 
element studies are presented in figures 3 to 5 for various values of mean line pressure. 
These experimental values were obtained by the method previously described in the sec- 
tion PROCEDURE. The analytical results (curves) are presented for comparison. The 
curves were obtained by use of equations (3) to (7). Curves for frequencies less than 
7 cps are shown as dashed lines because, for these frequencies, the distributed series 
disturbance -resistance may be significant in the computation of the attenuation and phase 
constants (appendix B). 

The attenuation constant a (in Np/ft) is presented in figure 4 as a function of the 
perturbation frequency. The agreement between the analytical and experimental results 
is considered good. The experimental data exhibit some scatter and slight variations 
with mean line pressure. Because of the high accuracy requirements, however, these 
effects can be attributed to limits of accuracy of the instrumentation. 

A comparison of the results for the porous -metal and lossy-tube shunt-element 
studies in figure 3 shows the effect of inertive impedance in the shunt elements. For the 
porous -metal element studies, a was essentially constant for frequencies greater than 
30 cps, while for the lossy-tube elements, a decreased with increasing frequency over 
this frequency range. This behavior occurred because for the porous -metal elements 
G d (eq. (5)) and were independent of frequency; whereas for the lossy-tube elements, 
the inertance had a greater effect on decreasing the value of G D (eq. (6)) than increasing 
K d (eq. (7)). 

Equation (3) shows that, for a fixed frequency and an essentially constant value of 
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Experimental 


(a) Porous-metal shunt element. Mean pressure, 75 pounds per square Inch gage (b) Porous-metal shunt element. Mean pressure, 155 pounds per square inch 


(5. 17x10 s N/sq m gage). 


gage (1.07X10 6 N/sq mgage). 





















y 




Tp- 

£> 










9 

/6 








~cr~ 


O 




/ 










o 5 





r 

ai 





'OQOC 

n: — ^ 

1 

— i- 



0 10 20 30 40 50 60 70 80 90 0 10 20 30 40 50 60 70 
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(c) Lossy-tube shunt element. Mean pressure, 97 pounds per square inch gage 
(6.69x10 s N/sq m gage). 


gage (9.51xl(r N/sq m gage). 

Figure 3. - Variation of attenuation constant with perturbation frequency. 


(d) Lossy-tube shunt element. Mean pressure, 138 pounds per square inch 
gage (9.51x10 s N/sq m gage). 



K d , a decreases as decreases. If die radius of each lossy tube were sufficiently 
small (inequality, (eq. (Cl)), a would not decrease with increasing frequency for the 
lossy-tube shunt-element case over the frequency range of consideration. 

In order to interpret the results obtained for a more easily, the attenuation con- 
stant is presented in figure 4 in terms of e _1 ^ 0a . This term is the ratio of downstream 
to upstream pressure for a wave traveling in the downstream direction for two points 
along the line that are separated by 100 feet (30. 5 m). For the same test line with no 
shunt elements (ref. 1), the attenuation is negligible. The results of figure 4 indicate 
that (except for the lower frequencies) the attenuation was approximately 55 percent for a 
100-foot (30. 5 m) line for the porous-metal shunt-element studies, and was an average of 
approximately 65 percent for a 100-foot (30. 5 m) line for the lossy -tube shunt-element 
studies. There is good agreement between the design (analytical) and experimental 
values. Therefore, it has been experimentally shown that sinusoidal perturbations in a 
fluid flowing in a long line can be attenuated to a design value. High values of attenuation 
may be realized by increasing the value of die effective distributed leakage Gp (eq. (3)). 
Equations (5) and (6) show that the value of °D can be indreased by (1) decreasing the 
lumped resistance of each shunt assembly R g (or R T ) (see appendix C), (2) increasing 
the number of shunt assemblies per unit length of main line n/£^, or (3) decreasing the 
inertance of each shunt assembly L T (lossy-tube shunt-element case). 

Analytical and experimental values of the phase constant 0 are presented in figure 5 
in terms of die ratio of phase constant to perturbation frequency as a function of fre- 
quency. The analytical and experimental results agree very well. 

As in reference 4, a comparison of the experimental values of the complex dimen- 
sionless downstream admittance z/Z ^ (eqs. (1) and (2)) with that calculated from the 
analogous electrical transmission line equations (ref. 5) and the resistance of the down- 
stream orifice showed no serious discrepancies. (Magnitudes agreed within ±10 percent, 
and angles agreed within ±10°. ) This agreement confirms the implied assumption that 
the data analysis gave the roots of interest. 

No strong conclusion was reached as to which of the two types of tested shunt ele- 
ments is more advantageous. During the determination of the leakage of the porous - 
metal shunt elements, it was noted that the elements could be partly clogged by dirt par- 
ticles rather easily. This fact may limit the practical application of porous-metal shunt 
elements. From the analysis of the lossy-tube shunt element, however, (eq. (3) and 
appendix C) the attenuation constant was found to be very sensitive to the diameter of the 
lossy tubes. A 1 -percent increase in the diameter of each lossy tube would cause about 
a 4 percent increase in the attenuation constant. Therefore, the lossy tubes must be con- 
structed carefully in order to obtain the design value of the attenuation constant. 
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Pressure ratio, 





Perturbation frequency, f, cps (Hz) 


(d) Lossy-tube shunt element. Mean pressure, 138 pounds per square inch gage 
(9.51X10 5 N/sq m gage). 

Figure 4. - Effect of frequency on attenuation, expressed as pressure amplitude ratio e~ 100a , (where 
a is attenuation constant. 
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(a) Porous-metal shunt element. Mean pressure, 75 pounds per square inch gage (b) Porous-metal shunt elements. Mean pressure, 155 pounds per square 

■ . {5. 17xl(r N/sq m gage). inch gage (1. 07xl(r N/sq m gage). 
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Perturbation frequency, f, cps (Hz) 


(c) Lossy-tube shunt element. Mean pressure, 97 pounds per square inch gage (d) Lossy-tube shunt element. Mean pressure, 138 pounds per square inch 

(6. 69xKP N/sq m gage). gage (9. 51xl(r N/sq m gage). 

Figure 5. - Effect of frequency on ratio of phase constant to perturbation frequency. 




SUMMARY OF RESULTS 


The results of this study show that (1) shunt assemblies closely spaced relative to 
sinusoidal disturbance wavelength can effectively attenuate to design values sinusoidal 
perturbations at a given frequency superimposed on the mean flow of a fluid in a long 
line, and that (2) the analogous electrical transmission line theory adequately predicts 
the attenuation constant when series disturbance -resistance is negligible. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 27, 1966, 

129-02-04-06-22. 



APPENDIX A 


SYMBOLS 


o 

A cross-sectional area of a line, ft 

(m 2 ) 

a isentropic (R^ = 0) sonic speed 
in the line with zero leakage, 

ft/sec (m/sec) 

Cp distributed shunt compliance of 
the main line not including the 
reactance (negligible for porous - 
metal shunt elements) due to 
shunt assemblies, ft 4 /lb (m 4 /N) 

d length of each lossy tube, ft (m) 

f frequency, cps (hertz) 

Gjj effective distributed leakage of 

the line, ft 4 /lb-sec (m 4 /N-sec) 

h length of line, ft (m) 

^ defined by eq. (C20) 

I defined by eq. (C21) 

> ^ 

K lumped compliance, ft^/lb (m^/N) 

K D effective distributed compliance of 
the line (C^ for porous -metal 
shunt element case), ft 4 /lb 
(m 4 /N) 

k constant defined by eq. (B7) 

£ tube length from any point to 
station C (fig. 1(b)), ft (m) 

2,5 

L lumped inertance, lb-sec /ft 
(N-sec 2 /m 5 ) 

L D distributed series inertance of the 
line, lb-sec 2 /ft 6 (N-sec 2 /m 6 ) 


L_ distributed specific inertance of 

T 2/4 

each lossy tube, lb-sec /ft 

(N-sec 2 /m 4 ) 

M constant defined by eq. (C5) 

n number of shunt assemblies (33) 

P complex amplitude of perturbation 

pressure, lb/ft 2 (N/m 2 ) 

AP pressure amplitude drop across 
length of lossy tube (complex), 
lb/ft 2 (N/m 2 ) 

p pressure, lb/ft 2 (N/m 2 ) 

Ap pressure drop between stations A 

and C (fig. 1(b)), lb/ft 2 (N/m 2 ) 

R lumped disturbance -resistance, 

lb-sec/ft^ (N-sec/m®) 

distributed series disturbance - 
resistance of the line, 
lb-sec/ft® (N- sec/m®) 

R t distributed specific disturbance - 

resistance of each lossy tube, 
lb-sec/ft 4 (N-sec/m 4 ) 

r inner radius of each lossy tube, 

ft (m) 

U complex amplitude of sinusoidal- 

perturbation velocity in axial 
direction of each lossy tube aver- 
aged over a cross-sectional area 
of each tube, ft/sec (m/sec) 
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V complex amplitude of sinusoidal- 

perturbation volume flow rate 

in axial direction for each 
2 3 

lossy tube, Uirr , ft /sec 
(m^/sec) 

v volume flow rate in the line, 

ftVsec (mVsec) 

Yp effective distributed shunt admit- 

tance of the line (complex), 
ft 4 /lb-sec (m 4 /N-sec) 

Z lumped impedance, lb-sec/ft^ 

(N-sec/m^) 

Zq impedance at station C (fig. 1(b)) 

looking downstream (complex), 
(approximately real for this 
experiment), lb-sec/ft^ 
(N-sec/m^) 

Z~ distributed series impedance of 

JJ 0 

line (complex), lb-sec/ft 
(N -sec/m®) 

z characteristic impedance of the 

line (complex), Ib-sec/ft^ 
(N-sec/m 5 ) 

a attenuation constant, neper s/ft 

(Np/m) 

j3 phase constant, rad/ft (rad/m) 

y ratio of specific heats, dimen- 

sionless 

6 defined by eq. (B14) 


p 

\i absolute viscosity, lb-sec/ft 

2 

(N-sec/m ) 

£ coordinate in axial direction of 

each lossy tube, ft (m) 

3 3 

p density of liquid, slugs/ft (kg/m ) 

a maximum radius of tubes in a 

series -parallel combination 
equivalent to a typical porous - 
metal shunt element, ft (m) 

cl> frequency, 2zrf , rad/sec 

0 volume of gas in shunt assembly 

3 3 

accumulator, ft (m ) 

Subscripts: 

A, B, C reference stations located as 
shown in fig. 1(b) 

g shunt-accumulator condition 

when charged 

1 case of = 0 

m mean value 

s porous -metal shunt assembly 

T lossy-tube shunt element 

0 zeroth order 

1 first order 


14 


■ I III 


■ nil i 


II II II I I 



APPENDIX B 


DERIVATION OF ANALYTICAL ATTENUATION AND PHASE CONSTANTS 

Attenuation of pressure disturbances due to factors such as heat exchange, absorp- 
tion of acoustic energy by the main tube wall, molecular exchanges of energy, and im- 
purities (ref. 6) are neglected in the following derivation of the attenuation and phase 
constants . 


DERIVATION OF EQUATIONS (3) AND (4) OF THE REPORT 

By analogy with electrical transmission line theory (ref. 6), the propagation constant 
is given by (ref. 5) 


a + j/3 = V z D y d CB1) 

where the effective distributed shunt admittance of the line is 

Y D = g d + 3“ K d ( B2 ) 

and, for the case of zero distributed series disturbance -resistance (R-q = 0), the distrib- 
uted series impedance is 


Z D = jwL D (B3) 

Substituting equations (B2) and (B3) into equation (Bl), squaring the result, and 
equating the real and imaginary parts yield 

/3 2 - a\ = co 2 L d K d (B4) 

and 

2/Sor ^ = wLpGjj (B5) 

o 

Solving equation (B4) for p and taking the positive square root of the result gives equa- 
tion (4) of the text. 

In order to find the expression for o^, equation (B5) is solved for j8 and the result 
is substituted into equation (B4). This procedure gives, after the terms are rearranged, 
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(B6) 


“i 4+ " 2l d k d“i -(j l d g d) 2 = ° 

2 

Solving this quadratic equation for o', and taking the positive square root of the result 
and rearranging terms give equation (3) of the text for the attenuation constant a 


SUBSTANTIATION FOR HIGHER FREQUENCIES OF ASSUMPTION THAT R Q = 0 


The equation for the distributed series -disturbance resistance of the line is 
obtained from the equation for the pressure drop in a line for turbulent flow as follows: 

AP = kv 2 (B7) 


It is assumed that 


r d ~ 


1 d(AP) 
1 A dv 


(B8) 


From equations (B7) and (B8) 


R = 2(AP) 
D v£ A 


(B9) 


This equation is often used for small disturbances superimposed on turbulent mean 
flow instead of a rigorously derived equation. Another equation that might be used to 
evaluate Rq for turbulent mean flow is the one for laminar disturbance flow (R„/A, 
eq. (C19)). Of these two equations, however, equation (B9) gave the greatest values of 
Rq for this experiment. Equation (B9) is, therefore, used as the worst case equation. 

For the nonzero values of distributed series line, disturbance -resistance, equa- 
tions (Bl) and (B2) are valid, but the following equation for is the more general 
equation instead of equation (B3): 

Z D = Rq + (BIO) 

o 

Following the procedure that led to of for the case of R^ = 0 presented previously 
gives 
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2 

The effect of Rq for the porous -metal shunt element case is now examined. The 
maximum value of Rp for the runs of this case (obtained from eq. (B9) with use of mean 

state data) is 9. 76xl0 2 (5. 41x10® N -sec/m®) pound -second per foot®. The values of 

G D’ L D’ and K D are S iven 111 a PPendix C. 

Using the numerical values of these line parameters and a frequency of 7 cps or 
greater gives 


r d g d " w2l d k d + wL D )l 


“ 2 4 


+ Gp { 1 + 


4 

co 2 L 2 


, 1/2 


(Bll) 



2 , 2 

a, L d 


< 0. 004 


(B12) 


Neglecting this term in equation (Bll), taking the square root of the result, and rear- 
ranging terms give 


a 




(B13) 


where the remaining correction term for distributed series disturbance -resistance is 


6 = 


Hd G D 

60 Lq ^> 2 K 2 + G 2 - o>K, 




D ) 


26f 


(B14) 


(See eq. (3) of the text. ) Evaluating 6 (always positive) for test frequencies of 7 cps or 
greater gives 


6 < 0. 14 (B15) 

From inequality (B15) and a binomial expansion of the factor (1 + 6)1/ 2 in equation (B13), 
it follows that 


a 




t i 

w k d + 


D 



(B16) 
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The bracketed term raised to the one -half power is the exact attenuation constant for zero 
distributed series disturbance -resistance R D (eq. (3)). The magnitude of the correction 
term dj in equation (B16) for the nonzero value of is less than 7 percent for fre- 
quenciesequal to or greater than 7 cps over the frequency range of the experiment. This 
error is within the limits of accuracy of the experiment. Therefore, Rj-j can be ne- 
glected for these frequencies. 

In computing the phase constant 0 for tested frequencies of 7 cps or higher, the 
magnitude of the maximum correction factor for the nonzero value of R^ can be shown 
to be negligible by a similar procedure. 

This method applied to the lossy-tube shunt-element case showed that R^ can also 
be neglected for tested frequencies of 7 cps or greater. 



APPENDIX C 


DISTRIBUTED LINE IMPEDANCE DUE TO SHUNT ASSEMBLIES 

The acoustic line used in this study is treated as an analogous electrical transmis- 
sion line, as is generally done. The acoustic analogous of electrical voltage, current, 
inductance, and capacitance are pressure, volume flow rate, inertance, and compliance, 
respectively, (ref. 6). 

Before the two types of shunt assemblies are analyzed separately, it is shown that 
each can be treated as distributed along the main line. This distribution is established 
by proving that there is at least one shunt assembly per quarter wavelength (ref. 6) as 
follows. For a perturbation frequency of 82 cps (the maximum experimental value) and 
an approximate wave propagation velocity of 3920 feet per second (1195 m/sec) (isentropic 
zero -leakage value of acoustic velocity), the minimum sinusoidal disturbance quarter 
wave-length for the experimental study was 12 feet (3. 66 m). The shunt assemblies were 
placed 2 feet (6. 1x10 m) apart. Because the minimum number of shunt assemblies per 
quarter wavelength over the frequency range of the experiment was at least 6, the shunt 
assemblies can be regarded as distributed. 


POROUS-METAL SHUNT ASSEMBLIES 
Impedance of Porous-Metal Shunt Element 

The impedance of a porous-metal shunt element is shown to be approximately resis- 
tive and frequency independent. Laminar disturbance flow through a porous-metal shunt 
element is assumed. As a conceptual aid, each shunt element is treated as a series- 
parallel combination of tubes of maximum radius a. Then the impedance of a shunt 
element is approximately resistive and frequency independent provided that (ref. 6) 

<J<+P± (Cl) 

▼ pO) 

For an absolute viscosity ju of 1. 57x10 ~ 5 pound -second per foot 2 (7. 52X10 -4 N - sec/m 2 ) 
(JP-4 fuel at 79° F), a frequency co of 80x2 u radians per second, and a density p of 1. 50 
slugs per foot 2 (7. 73x1 0 2 kg/m 2 ) (ref. 1), inequality (Cl) becomes (worst case value) 

ct < 1.44X10 -4 ft= 0. 438X10 -4 m (C2) 
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The manufacturer states that each porous -metal shunt element used filters out parti- 
cles with a maximum dimension greater than 4. 17x10"^ foot (1. 27x10"^ m). The right 
side of inequality (C2) is 6. 9 times this amount. If it is now assumed that the effective 
maximum tube radius of the pores in each shunt element a is within 6. 9 times the 
largest dimension of the filtered particles, it follows that the condition set by inequality 
(Cl) is satisfied. Therefore, the impedance of each porous-metal shunt element is ap- 
proximately resistive and frequency independent. 

For this reason, the value of the resistance R g can be determined from curves of 

static -pressure drop plotted against flow rate (nearly linear). Using this method yields 

a value of R of 4. 05x10^ pound -second per foot'* (6. 85x10*® N-sec/m^). 
s 


Comparison of Shunt-Element Resistance with Shunt-Assembly Impedance 

The impedance of each shunt assembly (fig. 2(a)) is 

Z s - R s * 4 L s - <P*> 

where each assembly is treated as a lumped parameter series acoustic circuit consisting 

of the resistance of the porous-metal shunt element R , the compliance of the accumu- 

s 

lator K , and the inertance of the mass of fluid L . It will be shown that the value of 

b S 

Z g is approximately that of R g over the frequency range of the experiment. 


Compliance 

The compliance of the accumulator (from ref. 6) is 

K = — ^ (C4) 

S dp 

For a perfect gas and an adiabatic process 

pf2 r = M (C5) 


Equations (C4) and (C5) give 
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(C6) 


1 Si 

y 


Taking the gas volume Si and the absolute pressure of the gas in the accumulator p 
as the mean volume Si m and absolute mean test line pressure p m , respectively, gives 

Si 

K = -2- (C7) 

^m 

The accumulator was brought from its initial charged state to its mean test state iso- 
thermally. This procedure yields 


= 7 “ g < C8 > 

p m 

Substituting equation (C8) into equation (C7) gives 

Si p 

K = -JL&. (C9) 

s 2 
^m 

The minimum value of shunt-accumulator compliance calculated from equation (C9) for 
the runs was 


K g 1.22X10' 7 ft 5 /lb = 7.21X10 11 m 5 /N (CIO) 

s 

The minimum frequency used in the experiment was 0. 5 cps. Use of this value 
yields the maximum magnitude of impedance due to the compliance of the accumulator as 

— < 2. 61X10 6 lb-sec/ft 5 = 4. 41X10 9 N-sec/m 5 (Cll) 

u>K s 

The ratio of this impedance to the resistance of a typical porous-metal shunt element is 
-2 

less than 6. 5x10 over the frequency range studied; for frequencies of 7 cps or greater, 
its maximum value is less than 4. 7X10"^. Therefore, the compliance of the shunt accu- 
mulator K can be neglected, 
s 
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Inertance 


From reference 6, the inertance due to nonviscous fluid of density p in a line of 
length h and cross-sectional area A is 

L = p — (C12) 

A 

Adding the maximum inertance of each section of a single shunt assembly in series 
(excluding that in the porous -metal disk already shown to be negligible) and using 80 cps 
indicated that the magnitude of the impedance due to the inertance of the shunt assembly 
is 


wL < 6. 9X10 5 lb-sec/ft 5 = 1. 17X10 9 N-sec/m 5 (C13) 

s 

It follows that the ratio of the inertive impedance to the resistance of a shunt assembly is 
2 

less than 1. 8x10 over the frequency range tested. Therefore, the inertance of each 

shunt assembly L can be neglected, 
s 


Distributed Line impedance 


The impedance of each porous -metal shunt assembly is essentially resistive. 
Therefore, from the values given herein 



n 



1.24X10" 8 ft 4 /lb-sec * 2.40x10 11 m 4 /N-sec 


(C14) 


L n = — (see eq. (C12)) = 3. 62X10 2 lb-sec 2 /ft 6 = 2. OlxlO 6 N-sec 2 /m 6 (C15) 
u A 

and the effective distributed compliance of the line is the same as the distributed 
compliance for the line with no shunt assemblies C^. From reference 6 


C D ~ 


L D a 


1. 80x10 10 ft 4 / lb = 3. 49X10 ” 13 m 4 /N 


(Cl 6) 


LOSSY-TUBE SHUNT ASSEMBLIES 


Impedance of Lossy-Tube Shunt Elements 

In order to calculate the impedance of a lossy-tube shunt element, the disturbance 
flow through a lossy tube is assumed to be laminar. For laminar disturbance flow, ref- 
erence 4 gives 


_££=RU+jTLU (C17) 

di 7 7 

where (1) U is the complex amplitude of the sinusoidal-perturbation velocity in the axial 
direction averaged over a cross-sectional area of a tube, (2) the distributed specific in- 
ertance of one tube is 



(Cl 8) 


(3) the distributed specific disturbance -resistance of one tube is 



(C19) 


and (4) I x and I y are expressed in terms of the Bessel-Kelvin functions of the first 
kind as, respectively, 
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The complex amplitude of volume flow rate at a point along the length of a tube is 


V = U?rr 2 

Combining equations (C17) and (C22) gives 



(C22) 


(C23) 


The ratio of tube length d to wavelength was less than 0. 003 for the frequencies of 
the experiment. For this reason the right side of equation (C23) is taken as independent 
of (This is equivalent to taking the distributed shunt compliance of each lossy tube 
as zero. ) The difference in pressure amplitude across the length of a tube at a given 
frequency is then 


'Rd L d 
AP =( — 2_+ jo -1—1 V 


(C24) 


, jrr 


nr . 


Each shunt element contained 21 tubes in parallel. From equation (C24), the lumped in- 
ertance of a shunt element is 


L r d 

L t = — — (C25) 

2177T 2 

the lumped resistance is 

Rd 

R T = — — (C26) 

1 2 
21j7T 

and the lumped impedance is 

Z T = R t + jo>L T (C27) 

The lumped inertance Lrj, and the lumped resistance R^, of a lossy -tube shunt el- 
ement can be evaluated for a given frequency by using tables of ber and bei functions 
(ref. 7). The values of the parameters used to obtain Lrj, and R,p were 


24 



p = 1 . 50 slugs/ft 3 = 7. 73x1 0 2 kg/m 3 
r = 3. 06X10 -4 ft = 0. 933x10 ” 4 m 

li = 1. 54X10 " 5 lb -sec/ft 2 (JP -4 fuel at 85° F) = 7. 38x10 “ 4 N -sec/m 2 
d = 0. 125 ft = 3. 81'XlO 2 m 

The lumped resistance and inertance were essentially constant over the frequency 
range of the experiment with the following values, respectively: 

R t = 2. 66X10 7 lb-sec/ft 5 = 4. 50X10 10 N-sec/m 5 


and 


L t = 4. 04X10 4 lb-sec 2 /ft 5 = 6. 84X10 7 N-sec 2 /m 5 

A calculation using the method given in the porous-metal shunt-element section of 
this appendix showed that the impedance of the rest of the lossy-tube shunt assembly 
(accumulator, tubing, etc. , fig. 2) was negligible compared with that of the lossy -tube 
shunt elements. 


Distributed-Line Impedance 

From equation (C27) the effective distributed shunt admittance including the shunt 
compliance of the main line is 


y D “ G D + J wK d 

where the effective distributed leakage is 


g d = — — 
D t.2 


A R“ + o> 2 L 2 


(C28) 


(C29) 


and the effective distributed shunt compliance of the line is 


K D “ C D 


n 

SL 


a o2 a „2 t 2 

A Hjp ^ w Lip 


(C30) 


The distributed inertance of the line is the same as that for the porous -metal shunt 
element case (eq. (C15)). The value of is given by equation (Cl 6). 
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